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Abstract 

We prove the Euler-Lagrange delta-differential equations for problems of the calculus of variations 
on arbitrary time scales with delta-integral functionals depending on higher-order delta derivatives. 
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1. Introduction 



In recent years numerous works have been dedicated to the calculus of variations on time scales 
and their generalizations — see 0, Il2 



H El, ED, 0, m 




and the references therein 



Most 

of them deal with delta or nabla derivatives of first-order 12, 3 only a few 

with higher-order derivatives nam. Depending on the type of the functional being considered, 
different time scale Euler-Lagrange type equations are obtained. For variational problems of first- 
order the Euler-Lagrange equations are valid for an arbitrary time scale T, while for the problems 
with higher-order delta (or nabla) derivatives they are only valid in a certain class of time scales, 
more precisely, the ones for which the forward (or backward) jump operator is a polynomial of 
degree one da Hi]. Here we consider variational problems involving Hilger derivatives of higher 
order, and prove a necessary optimality condition of the Euler-Lagrange type on an arbitrary time 
scale, i.e., without imposing any restriction to the jump operators. 



2. Preliminaries 

Here we recall some basic results and notation needed in the sequel. For the theory of time 
scales we refer the reader to 

A time scale T is an arbitrary nonempty closed subset of the real numbers R. The func- 
tions a : T— >T and p : T— >T are, respectively, the forward and backward jump operators: 
a(t) = inf {s £ T : s > (} with inf = supT (i.e., a(M) = M if T has a maximum M); p(t) = 
sup{s G T : s < i} with sup0 = inf T (i.e., p(m) = m if T has a minimum m). The symbol 
denotes the empty set. The graininess function on T is defined by p(t) := a(t) — t. For T = R 
one has a(t) = t = p(t) and p{t) = for any t G R. For T = Z one has cr(t) = t + 1, p(t) = t — 1, 
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and n(i) = 1 for every t £ Z. A point t £ T is called right-dense, right-scattered, left-dense, or 
left-scattered, if a(t) = t, ait) > t, pit) — t, or pit) < t, respectively. 

Let T = [a, b] n T with a < b and T a time scale. We dehne T K := T\ (/>(&), 6], and TT° := T, 
T K " := ^T K " J for n £ N. The following standard notation is used for a (and p): a°(t) = t, 

a n (t) = (o-oer"- 1 )^), n £ N. 

We say that a function / : T — ► R is delta- differentiable at £ £ T K if there is a number / A (<) 
such that for all £ > there exists a neighborhood U of t such that 

|/(er(t)) - /(s) - / A (t)(cr(t) - s)| < £|cr(i) - s\, for all s £ U. 

We call f A it) the delta-derivative of / at i. We note that if the number / A (i) exists then it is 



unique in T K (see 1% |l5[). In the special cases T = R and T = Z, / A reduces to the standard 
derivative f'{t) and the forward difference A/(i) = /(f + 1) — /(£), respectively. Whenever / A 
exists, the following formula holds: fit) = fit) +/i(i)/ A (t), where we abbreviate foa by f a . Let 

/ A ° = /. We dehne the rth-delta derivative of / : T K ' -> M, r £ N, to be the function (/ Ar_1 ) ; 

provided / A is delta differentiable on T K . 

A function / : T — > K is called rd-continuous if it is continuous at the right-dense points in T 
and its left-sided limits exist at all left-dense points in T. A function / : T — > K™ is rd-continuous 
if all its components are rd-continuous. The set of all rd-continuous functions is denoted by C r d- 
Similarly, C"! rd will denote the set of functions with delta derivatives up to order r belonging to 
C r d- A function / is of class / £ C r prd if f A is continuous for i = 0, . . . , r — 1, and / A exists and 
is rd-continuous for all, except possibly at finitely many t £ T K . 

A piecewise rd-continuous function / : T — > R possess an antiderivative F A = /, and in this 
case the delta integral is defined by f f(t)At — Fid) — F{c) for all c, d £ T. It satisfies 

/(r)Ar = mW/W. 

Jt 

b b 

If T = R, then J fit) At — J f(t)dt, where the integral on the right hand side is the usual Riemann 

a a 

b b-l 

integral; if T = Z and a < b, then / fit) At = £ f(k). 

a k—a 

3. Main results 

Consider the following higher-order problem of the calculus of variations up to order r, r > 1: 

fP r - 1 (b) 

Ciyi-)) = / Lit, yit), y A it), ...,y A it)) At — ► min, (1) 

J a 

subject to boundary conditions 

y(a) = y° a , v{p r -\b))=yl,-.- ,y Ar -\a)=y r -\ y A ^ {p r -\b)) = y\-\ (2) 

where T is a bounded time scale with a :— minT and b := maxT, L : [a, p r ib)]f x R r+1 — >• R is 
a given function, where we use the notation [c, d]y :— [c, d] n T, and y a , y\ £ R, % = 0, . . . , r — 1. 
The results of the paper are trivially generalized for functions y : [a, b]r — > R n , but for simplicity 
of presentation we restrict ourselves to the scalar case n = 1. 

A function y(-) £ C r prd is said to be admissible if it is satisfies condition (|2|). An admissible 
y(-) is a weak local minimizer for ©-(HJ if there exists 6 > such that Ciyi-)) ^ £(z/(')) f° r an y 
admissible y £ Cp rd with \\y — y|| r ,oo < ^, where 

r 

llylkoc := V 2/ A * 

II oo 

i=0 
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y A — y and ||y||oo : = su Pte[a.p>-(b)] T For simplicity of notation we introduce the operator 

[y] defined by [y](t) = (t, y{t) 1 y A (t), . . . , y Ar (t)) . Then, functional (fTJ) can be written as 

C{y{-)) = / L[y](t)At. 

We assume that (t,u±, . . . , u r+1 ) — > L(t, u\, . . . , w r +i) has continuous partial derivatives for 

all t € [a,p r (b)]j, i = 1, . . . ,r + l, and i — >• and t — > ^-\y]{t), i = 1, . . . ,r + l, are piecewise 

rd-continuous for all admissible functions y(-). 



3.1. The higher-order Euler- Lagrange equation 

We now prove the Euler-Lagrange equation for problem ([l])-©. 

Remark 1. In order for the problem to be nontrivial we require the time scale T to have at 
least 2r + 1 points. Indeed, if the time scale has only 2r points, then it can be written as 
T = {a,a(a),...,a 2r - 1 {a)} and 

/ L(t,y(t),y A (t),... 7 y Ar (t))At 

= ^ L(t,y(t),y A (t),...,y Ar (t))At = J2 f ^ L(t,y{t),y A (t), . . . ,y AT (t))At 

r-1 

= $> i+1 (a) - a\a))L(o\a),y{o\a)), y A (a»), . . . l2/ A »)))- (3) 



Having in mind the boundary conditions and the formula f A (t) = , we can conclude 

that the sum in ([3]) is constant for every admissible function y(-). 

Theorem 1. If y(-) is a weak local minimizer for the problem (U}-©, ^hen y(-) satisfies the 
Euler-Lagrange equation 



dL . ... 

—IvW) 



Oy 



(-iy 



dy 



i=0 

a(t) ( r a(r r ) 



a J a 



— Ar _ 2 _i [y](T r -i-i)AT r _i_i • • • Ar r _iAr r 



dy L 



— [y](T 1 )A Tl + ci 
oy 



Ar r _ 1 — (— l) r c r _i > Ar r — c r = 



/or some constants C\, . . . , c r and all t G [a, ( o r (&)]x. 

Proof. We first introduce some notation: yo(t) = y(t), Vi(t) = D A {t), 
u(t) = y A (t). Then problem ([T])-© takes the following form: 



(4) 



y r _i(t) = y A ~ (t) 



£[»(•)] = / £(*, 2/o (t), Vi (t), • • • , y r -i(t),u(t))At — > min, 

■/ a 

yf(t) = y i+1 (t), i = 0,...,r-2, 

y r A 1 (*) = «(*), 

^(a) = tf, ^ (p'-^ft)) = ?4, j = 0, . . . ,r - 1 . 
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With the notation x = (yo, yi, ■ ■ ■ , y r -i), our problem (HJ) — d2J) can be written as the optimal control 
problem 

L(t, x(t), u(t))At — > min, 



C[x(.)] = 



x A (t) = Ax(t) + Bu{t) . 



(5) 



where 



( 1 
1 



.4 



x(a) - x a 
x^ib)) - x b 








\ 



1 

/ 



B 



Note that assumption (Al) of 13, Theorem 9.4] holds: matrix / + fJ,(t)A is invertible, and the 
matrix r Vip(x(a), x(p r ~ 1 (b)) has full rank. Therefore, if (x(-),u(-)) is a weak local minimum for 
([5]), then there exists a constant A and a function p : [a, p r ~ 1 (b)]f -> M r , p G C^ rd , such that 
(A, p(-)) ^ and the following conditions hold: 



-p A (t) = A T p°(t) + \ 



dL 



— (t,x(t),u(t)) 



B 1 p°(t) + \—(t,x(t),u(t))=0 

for all t G [a, /9' (6)]t. Consequently, if y(-) is a weak local minimizer for ([T|)— (J2j) , then 

Fit 

tf-i(*) = -%M(t) 
holds for all t G [a, p r (b)]j, where p^_ 1 {t) is defined recursively by 



PoW = 
Pf(*) 

with Ci, i = 0, . 



(*) = -/ A— MCnJAn-d, 

Jo oy 



dL, 



Ar»+i— Ci_i, i = l,...,r-l. 
r — 1, constants. From ([T])-© we obtain that equation 



(6) 
(7) 

(8) 
(9) 



[y](r r )Ar,. 



r-3 



(-1)' 



+E(- 1 >' 

■i=0 
r(t) f rCT (r r ) 



<r(t) r <r(r r ) 



a J a 



A- [y](r r _i_i)Ar r _i_i • • ■ Ar r _iAr r 



er(r 2 ) ^ 

A o— [j/]( T l)Ari + Cl 



Ar r _i - (-l) r V-l } Ar r - c r = 



(10) 



holds for all t G [a, p r {b)]r. We show next that A ^ 0. First observe that if / G and f a (t) = 
for all t G [a, &]f , then f(t) = for all t G [c(a), 6]t- Suppose, contrary to our claim, that A = in 
equation ([5]) and (|7|). Then, we can write the system of equations 



=0, 

P?(t) =~PtM i = l,-..,r-l 
P?_i(t) =0, 



(11) 
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for all t G [a,p r (b)]j. From the last equation we have p r -i(t) — 0, Mt G [a(a), p r ~ 1 (b)]f. This 
implies that p A _ 1 (t) = 0, Vi G [c(a), p r (&)]T, and consequently Pr-2^) = 0, Vi G [cr(a), j o r (&)]T- 
Therefore, p r -2(t) = 0, Vi G [c 2 (a), j o r_1 (6)]T- Repeating this procedure we have pi(£) = for 
all t G [a r - 1 (a),p r - 1 (fi)]r. Hence, = = = -p A (t)n(t) - p (t) = -p (t) for all 

i G [cr r_1 (a), p r (&)]T- Note that the first equation of (jlll) implies po(t) = c for some constant 
c and all i G [a, //~ 1 (6)]x. Since the time scale has at least 2r + 1 points (see Remark [lj, the 
set t G [<7 r_1 (a), p r ^ 1 {b)]f is nonempty and we conclude that po(t) — for all t G [a, j o r ~ 1 (&)]i'. 
Substituting this into the second equation we get p A (t) = d for some constant d and all t G 
[a, p r ~ 1 (b)]j. Having in mind that pi(to) = for some to G [a, p r ~ 1 (b)]j we obtain pi(t) = for 
all t G [a, // -1 (&)]t- Repeating this procedure we conclude that Pi(t) = 0, i = 1, . . . , r — 1, for all 
t G [a, ( o r ~ 1 (fo)]T- This contradicts the fact that (A, £>(•)) ^ 0. Hence, equation (|TU)) can be divided 
by A and ((U) is proved. □ 

3.2. Corollaries 

For illustrating purposes we consider now the two simplest situations, i.e., r — 1 and r = 2. 
Corollary 1 (cf. (5, [l6j|). Ify(-) is a weak local minimizer for the problem 

,b 

C{y{-))= L(t,y(t),y A (t))At ^ min 



subject to boundary conditions y{a) = y a and y(b) = yb, then y(-) satisfies the Euler- Lagrange 
equation 

H (t,y(t),y A (t)) = ( W |(^W,«/ A (r))Ar + Cl 
for some constant c\ and all t G [a, 

Corollary 2 (cf. [ljj H^]). If y(-) is a weak local minimizer for the problem 

CP(b) 



£(*/(•)) = / L(t,y(i),2/ A (t),y AA )Ai — > min 

J a 

subject to boundary conditions y(a) = y®, y(p(b)) = yb, y A (a) — y\, and y A {p{b)) — y\, then y(- 
satisfies the Euler- Lagrange equation 

)L (t,y(t),y A (t),y AA (t))- / ^(T2,y(T2hy A (T 2 ),y AA (T 2 ))AT 2 



% Mnw,s w ' y y " J a 9y A 

,<r(t) 



CT ( 1 "2) ar 

^-(r 1: y(T 1 ),y A (r 1 ) 7 y AA (r 1 ))Ar 1 + c 1 



Ar 2 - c 2 = 



/or some constants c\ and c 2 and all t G [a,p(6)]^. 



3.3. An example 

Let T = [a, &] n /iZ, where hi :— {hz\z G Z}, h > 0. Then for any / G Cp rd we have 



/ At 



/ AV \ iG{0,...,i-l}, (12) 



j — i integrals 

where / A CT (t) stands for f A (a^~ % (t)). We will show this by induction. For j = 1 
( %(£)A£= /*/(C)A^+ [ t+h f(0A£= f f(Z)A£ + hf(t), 
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and then 



Then for j = k + 1 



f° W / (£)A£ = f(t) + hf A (t) = f a . Now assume that dT2j) is true for all j = 1, . . . , fc. 



A 



r(t) 



a a 



/ At 



/At + ft 



a(t) 



/At 



k-j-l—i integrals 



k+l—i 



k—i 



\ k-i 



/At 



/ 



r(t) 



/At 



V 



k — i 



J 



Delta differentiating r times both sides of equation (Q} and in view of (|T2l , we obtain the /i-Euler- 
Lagrange equation in delta differentiated form: 



Lfa (t, y, y A , . . . , y Ar ) + ^-l)'"^?' - ' (t, y, y A . . . , y A ") = 0. 



i=0 
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